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ABSTRACT 

Eight destinations in the Solar System have sufficient 
atmosphere for aeroentry, aeroassist, or 
aerobraking/aerocapture: Venus, Earth, Mars, Jupiter, 
Saturn, Uranus, and Neptune, plus Saturn’s moon Titan. 
Engineering-level atmospheric models for Earth, Mars, 
Titan, and Neptune have been developed for use in 
NASA’s systems analysis studies of aerocapture 
applications. Development has begun on a similar 
atmospheric model for Venus. An important capability 
of these models is simulation of quasi-random 
perturbations for Monte Carlo analyses in developing 
guidance, navigation and control algorithms, and for 
thermal systems design. Characteristics of these 
atmospheric models are compared, and example 
applications for aerocapture are presented. Recent Titan 
atmospheric model updates are discussed, in 
anticipation of applications for trajectory and 
atmospheric reconstruct of Huygens Probe entry at 
Titan. Recent and planned updates to the Mars 
atmospheric model, in support of future Mars 
aerocapture systems analysis studies, are also presented. 

1. INTRODUCTION 

Engineering-level atmospheric models have been 
developed, or are under development, for five of the 
eight possible Solar System destinations where 
aerocapture could be used. These include Global 
Reference Atmospheric Models (GRAMs) for Earth 
(GRAM-99) [1, 2], Mars (Mars-GRAM 2001) [3-6], 
Titan (Titan-GRAM) [7], Neptune (Neptune -GRAM) 
[8], and Venus-GRAM (under development). Physical 
characteristics of the various planetary atmospheres 
vary significantly. Likewise, significant variation is 
found in the amount of available data on which to base 
the respective engineering-level atmospheric models. 
The detailed characteristics of these models differ 
accordingly. 

Earth-GRAM is based on climatology assembled from 
extensive observations by balloon, aircraft, ground- 


based remote sensing, sounding rockets, and satellite 
remote sensing. Details are provided in the GRAM 
User’s Guide [1]. Mars-GRAM is based on 
climatologies of General Circulation Model (GCM) 
output, with details given in the Mars-GRAM User’s 
Guide [3]. Mars-GRAM has been validated [4-6] by 
comparisons against observations made by Mars Global 
Surveyor, and against output from a separate Mars 
GCM. In contrast, data used to build Titan-GRAM and 
Neptune-GRAM are more limited, deriving primarily 
from Voyager observations and limited ground-based 
stellar occultation measurements. Titan-GRAM is based 
on data summarized in [9], while Neptune-GRAM was 
built from summaries of data contained in [10]. For 
Venus, a substantial amount of data has been collected 
from orbiter and entry probe observations. These have 
been summarized in the Venus International Reference 
Atmosphere (VIRA) [11], which forms the basis for 
Venus-GRAM (under development). 

Fig. 1 shows the wide variety of temperature profiles 
encountered among the planets and Titan. For Earth, 
Venus, Mars, and Titan, height is measured from a 
reference surface (mean sea level on Earth). On 
Neptune, height is measured above the level at which 
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Fig. 1. Comparison of temperature profiles among the 
planets and Titan. 
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pressure is one bar (Earth normal sea-level pressure). 
All of the planets exhibit a troposphere region, where 
temperature decreases with altitude, indicative of heat 
flow upward from the surface (on average). All of the 
planets exhibit a thermosphere region, where (on 
average) temperature increases with altitude, because of 
absorption of heat flux from the Sun as it penetrates into 
the atmosphere. All of the planets have stratospheres, 
where temperature decrease above the surface 
diminishes, and remains relatively constant until the 
base of the thermosphere (Earth being the exception to 
this, where the presence of ozone and resultant 
atmospheric heating produces a local temperature 
maximum in Earth’s stratosphere-mesosphere region). 

For interest in aerocapture or aerobraking, atmospheric 
density is the most important parameter. Fig. 2 
compares density profiles on the planets and Titan. 


Atmospheric Density Comparison 
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Fig. 2. Comparison of density profiles among the 
planets and Titan. 

Vertical dashed lines in Fig. 2 indicate typical density 
values at which aerocapture or aerobraking operations 
would occur. Intersections of the aerocapture dashed 
line with various density curves shows that aerocapture 
would occur at a wide range of altitudes at the various 
destinations, varying from about 50 km at Mars to about 
300 km at Titan. Aerobraking at Earth, Mars, and Venus 
would take place near, and just above, the 100 km level. 
At Neptune and Titan, aerobraking would be 
implemented near 550 km and 750 km, respectively. 

Fig. 2 shows that density decreases fairly rapidly with 
altitude for the terrestrial planets (Venus, Earth, Mars), 
while it decreases rather slowly for Neptune and Titan. 
This effect is explained by differences in density scale 
height, //, for the various planets and Titan. Density 
decreases rapidly with altitude if H is small, while it 


decreases slowly if H is large. H is proportional to 
pressure scale height [ R T / ( M g) ]. For the terrestrial 
planets, molecular mass M is large (M « 29-44), so H is 
small. On Neptune, H is large because M is small for 
Neptune’s hydrogen-helium atmosphere (M « 2). For 
Titan, H is large despite the high molecular mass of its 
atmosphere (M « 29), because its gravity is low. 

2. BASIS FOR THE ATMOSPHERIC MODELS 

In Earth-GRAM, Mars-GRAM, and Venus-GRAM, 
input values for date, time, latitude, longitude, etc. are 
used to calculate planetary position and solar position. 
In this manner, effects of latitude variation and seasonal 
and time-of-day variations can be computed explicitly. 
A simplified approach is adopted in Titan-GRAM and 
Neptune-GRAM, whereby these effects (as well as 
effects of relatively large measurement uncertainties for 
these planets) are represented within a prescribed 
envelope of minimum-average-maximum density versus 
altitude. Fig. 3 shows this envelope for Titan. 


Huygens (Yelle) Titan Atmosphere Model Density 
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Fig. 3. Minimum, average, and maximum density 
profiles for Titan [9]. 

Engineering atmospheric model data developed for the 
Huygens entry probe [9] are used to define the Titan 
envelope. For Neptune, data from [10] are employed to 
generate a comparable minimum-maximum envelope, 
as shown in Fig. 4. A single model input parameter, 
Fminmax, allows the user of Titan-GRAM or Neptune- 
GRAM to select where within the min-max envelope a 
particular simulation will fall. Fminmax = -1, 0, or 1 
selects minimum, average, or maximum conditions, 
respectively, with intermediate values determined by 
interpolation; i.e., Fminmax between 0 and 1 produces 
values between average and maximum. Effects such as 
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variation with latitude along a given trajectory path can 
be computed using the appropriate representation of 
Fminmax variation with latitude. 


Neptune Atmosphere Model Density Profiles 
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Fig. 4. Minimum, average, and maximum density 
profiles for Neptune from data in [10]. 

Since drag is proportional to density, density is the most 
important atmospheric parameter for aerocapture. Next 
most important is height variation of density (as 
characterized by density scale height). Density scale 
height is important in determining aerocapture corridor 
width, or entry angle range that allows the vehicle to 
achieve capture orbit without “skipping out” or 
“burning in”. As discussed above, small density scale 
height means rapid change of density with altitude, 
which results in low corridor width. Large density scale 
height implies slow density change with altitude, and 
large corridor width. 

Fig. 5 compares height profiles of density scale height 
among the planets and Titan. Aerocapture altitude (c.f. 
discussion of Fig. 2) is indicated by letter A in Fig. 5. 
This figure shows low density scale heights (4-8 km) at 
aerocapture altitudes for the terrestrial planets. Larger 
scale heights (« 30-50 km) occur at aerocapture 
altitudes on Neptune and Titan. 

3. TITAN-GRAM GCM OPTION 

An option has recently been added for using Titan 
General Circulation Model (GCM) data as input for 
Titan-GRAM. The Titan GCM data used are from 
graphs in [12]. Upper altitudes for the Titan GCM 
option are computed using a parameterized fit to Titan 
exospheric temperatures, taken from graphs in [13]. 


Density Scale Height Comparison 



Fig. 5. Comparison of atmospheric density scale height 
among the planets and Titan. 

Fig. 6 shows a height-latitude cross section of density, 
expressed as percent deviation from the mean, for 
Voyager encounter date November 12, 1980 

(planetocentric longitude of Sun Ls = 8.8°), 00:00 
GMT, longitude zero, local solar time 0.7 Titan hours. 
Fig. 7 compares vertical density profiles at latitude zero, 
local solar time 1 hour and 13 hours on the Voyager 
encounter date, with the Huygens Yelle [9] minimum- 
maximum density envelope from Fig. 3. This figure 
shows that the Titan GCM results correspond fairly 
closely with Yelle maximum conditions up to about 300 
km altitude, and agree quite closely with Yelle average 
conditions (vertical line at 0 in Fig. 7) above about 500 
km. 

Titan-GRAM GCM Option Density, % from Yelle Mean 
Voyager Time: 11/12/1980 00:00 Lon=0 Ls=8.8 LST=0.7Hrs 
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Fig. 6. Density (percent deviation from mean) versus 
height and latitude, using Titan-GRAM GCM option. 
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Yelle Density Envelope versus GCM 



Fig. 7. Comparison of two selected Titan-GRAM 
density profiles (GCM option) with minimum- 
maximum envelope from Huygens Yelle model [9]. 


4. VENUS-GRAM DEVELOPMENT 

Based on the Venus International Reference 
Atmosphere (VIRA) [11], Venus-GRAM is being 
developed and applied in ongoing Venus aerocapture 
performance analyses. Fig. 8 shows a plot of density 
(percent deviation from the mean) versus height and 
latitude from Venus-GRAM. Conditions in Fig. 8 are 
for Ls = 90° and local solar time =12 Venus hours. 


Venus Density, Percent from Average 
Ls = 90 LST = 12 Hr 



Fig. 8. Example height-latitude density cross section 
from Venus-GRAM. 

Below about 100 km altitude on Venus, we find that 
temperature, density, and density scale height 
conditions are very uniform with both latitude and time 


of day. VIRA data below 100 km altitude vary only 
slightly with latitude and have no dependence on local 
solar time. Between 100 km and 150 km, VIRA data 
depend on local solar time, but not latitude. From 150 
km to its top at 250 km, VIRA depends on solar zenith 
angle, which is affected by both latitude and local solar 
time. 

5. NEW MARS-GRAM FEATURES 

During Mars Global Surveyor (MGS) aerobraking 
operations, large density variations were observed 
between successive periapsis passes [14]. These 
appeared to be longitude-fixed or terrain-fixed waves, 
usually dominated by wave-2 or wave-3 components 
(wave-n meaning that n wavelengths fit around a 360° 
longitude circle). During Mars Odyssey aerobraking, 
similar large-amplitude density variations were 
observed. However, during some periods, Odyssey- 
observed density variations appeared to be traveling 
waves whose phase speed relative to a fixed longitude 
seemed to remain constant for a matter of a few days. 
Mars-GRAM 2001 has an option to represent terrain- 
fixed waves of the type observed by MGS. Work is 
underway to develop a new version of Mars-GRAM 
that will (among other features) include the option to 
allow user input values for phase speed of traveling 
wave components, of the type observed by Mars 
Odyssey. 

Also during MGS and Odyssey aerobraking, it was 
observed that Mars-GRAM produced better 
correspondence with observed atmospheric density if 
the altitude scale of its input Mars Thermospheric 
General Circulation Model (MTGCM) data base was 
shifted (described as a “height offset”). New sets of 
Mars General Circulation Model (MGCM) and 
MTGCM data are being produced for use as input in the 
next Mars-GRAM update. These GCM model runs 
include better treatment of the matchup conditions (both 
mean conditions and upward wave fluxes) between the 
upper boundary of MGCM and lower boundary of 
MTGCM (at the 1.32 pbar level, near 80 km). A new 
non-local thermodynamic equilibrium (non-LTE) 
method for treating near-infrared heating and C0 2 15- 
micron cooling will also be employed in the MTGCM 
model runs. This methodology is based on a non-LTE 
model of Lopez-Valverde and Lopez-Puertas [15]. 
More realistic dynamics in both MGCM and MTGCM 
data sets is also anticipated from the use of latitude and 
seasonal variations of dust optical depth observed by 
MGS Thermal Emission Spectrometer (TES) in its 
mapping years 1 and 2. It is hoped that these new 
MGCM/MT GCM input data sets for Mars-GRAM will 
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significantly lessen the need for height offset, and 
significantly improve the correspondence with observed 
densities during Mars-GRAM use in support of 
aerobraking operations for Mars Reconnaissance 
Orbiter. 

6. NEW MARS-GRAM SLOPE-WIND FEATURE 

For potential applications in preliminary site screening 
for Mars landers, a new slope wind feature is being 
developed for Mars-GRAM. Slope winds are computed 
in Mars-GRAM from a diagnostic (algebraic) 
relationship based on [16]. This approach differs from 
mesoscale models, such as Mars Regional Atmospheric 
Model System (MRAMS) [17], and Mars Mesoscale 
Model version 5 (MMM5) [18], which use prognostic, 
full-physics solutions to the time- and space-dependent 
differential equations of motion. As such, slope winds 
in Mars-GRAM will be consistent with its 
“engineering-level” approach, and will be extremely 
easy and fast to evaluate, compared with mesoscale 
model solutions. Mars-GRAM slope winds are not 
being suggested as a replacement for more 
sophisticated, full-physics mesoscale models, but may 
have value, particularly for preliminary screening of 
large numbers of candidate landing sites for future Mars 
missions. 

Terrain slopes used in the slope wind model are 
computed from 0.5° _ 0.5° Mars Orbiter Laser Altimeter 
(MOLA) topography. Mars-GRAM slope winds will be 
added to winds from MGCM, which have a resolution 
of 7.5° _ 9° in latitude and longitude. The Mars-GRAM 
slope wind model will thus add significantly higher 
resolution information about possible near-surface 
winds than is provided by MGCM. 

Fig. 9 shows Mars-GRAM slope winds, evaluated at a 
level 2 km above local terrain height for the Gusev 
Crater area, at the date and time of Rover Spirit landing. 
If this wind field is valid, then Spirit would have 
experienced up to ~ 25 m/s winds “opposing” its entry 
into Gusev Crater near an altitude of 1-2 km above 
surface level. Spirit experienced significant turbulence 
or winds during its descent (Prasun Desai, private 
communication), causing it to fire its Transverse 
Impulse Rocket System to correct for off- vertical firing 
of its main retrorockets, and to reduce its lateral impact 
speed. 

Fig. 10 shows MOLA terrain heights in a portion of the 
eastern end of Valles Marineris, used in these 
preliminary tests of the slope wind model. 


Gusev Wind 2 km Above Surface 



Fig. 9. Slope wind vectors at Gusev Crater, 2 km above 
surface, for date and time of Spirit landing. 


Terrain Height, km 



East Longitude, degrees 

Fig. 10. Terrain Heights in portion of Valles Marineris 
region. 

Fig. 11 shows northward component of Mars-GRAM 
slope winds, evaluated at a level 1 km above local 
terrain height for the study area shown in Fig. 10. The 
season assumed is Ls = 0° (northern spring equinox) at 
local time 14 hours. Comparison of Fig. 11 with Fig. 10 
shows that the major pattern for slope winds at this time 
is northward and upward along the north wall of the 
valley and southward and upward along the south wall 
(i.e. upslope flow on both valley walls), a reasonable 
situation for early afternoon local time. These examples 
of test output from the new Mars-GRAM slope wind 
model may be compared with wind simulations from 
Mars mesoscale models, presented by Rafkin and 
Michaels [19] and Kass, et al. [20]. 
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Northward Slope Wind, m/s 



East Longitude, degrees 


Fig. 11. Northward slope winds at Ls = 0° and LST = 

14 hours, 1 km above terrain surface. 

7. PERTURBATION MODELS 

An important feature of all the GRAM atmospheric 
models is their ability to simulate “high frequency” 
perturbations in density and winds, due to such 
phenomena as turbulence and various kinds of 
atmospheric waves. As illustrated in Fig. 12, Earth- 
GRAM altitude, latitude, and monthly variations of 
perturbation standard deviations are based on a large 
climatology of observations. 


Density Perturbations 



Altitude, km 


Fig. 12. Height variation of density perturbation model 
standard deviations for Earth, Mars, Titan, and Neptune. 

For Titan-GRAM and Neptune-GRAM, perturbation 
standard deviations are computed from an analytical 
expression for gravity wave saturation conditions, 
explained more fully in [7]. As shown in Fig. 12, the 
resulting vertical profiles of standard deviations for 
Titan and Neptune are not dissimilar to Earth 


observations, when expressed as percent of mean 
density. For Mars-GRAM, a similar gravity wave 
saturation relation is used to estimate density 
perturbation standard deviations, except that effects of 
significant topographic variation on Mars are also taken 
into account. Up to about 75 km altitude, the Mars 
model density standard deviations are also fairly 
consistent with Earth observations. By about 100 km to 
130 km altitude, Mars model density standard 
deviations increase to about 20% to 35% of mean value, 
consistent with observed orbit-to-orbit density 
variations observed by Mars Global Surveyor and Mars 
Odyssey. 

A typical application of the Neptune-GRAM 
perturbation model is shown in Fig. 13. 


Sample Neptune-GRAM Monte-Carlo Density Output 



Perturbed Density (% from Neptune Average) 


Fig. 13. Sample Monte Carlo density perturbations 
from Neptune-GRAM, expressed as percent deviation 
from Neptune mean value. 

Neptune-GRAM was recently utilized in Neptune 
aerocapture systems analysis studies. The chosen 
aerocapture design reference mission included 
simulations which involved capture into a highly 
eccentric orbit, to allow the orbiter to periodically visit 
Triton for scientific observations. The ability to 
successfully aerocapture into such an eccentric orbit 
depends very significantly on details of Monte Carlo 
trajectory simulations, particularly on atmospheric 
density variations such as illustrated in Fig 13. For such 
an eccentric orbit, there is relatively little margin for 
error between a captured orbit and one which exceeds 
escape velocity upon atmospheric exit, a result which 
could ultimately lead to mission failure. Neptune- 
GRAM was used to define an aerocapture corridor 
width consistent with mission success. 
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8. CONCLUSIONS 

The engineering-level atmospheric models presented 
here are suitable for a wide range of mission design, 
systems analysis, and operations tasks. For orbiter 
missions, applications include analysis for aerocapture 
or aerobraking operations, analysis of station-keeping 
issues for science orbits, analysis of orbital lifetimes for 
end-of-mission planetary protection orbits, and 
atmospheric entry issues for accidental break-up and 
burn-up scenarios. For lander missions to Venus, Mars 
and Titan, and for Earth-return, applications include 
analysis for entry, descent and landing (EDL), and 
guidance, navigation and control analysis for precision 
landing and hazard avoidance. Perturbation simulation 
capabilities of these models make them especially 
useful in Monte Carlo analyses for design and testing of 
guidance, navigation, and control algorithms, and for 
heat loads analysis of thermal protection systems. 
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